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l. INTRODUCTION 

A typical wind model used for ascent vehicle wind biasing consists 
of the monthly mean wind at each altitude. Such a model does not contain the 
small scale perturbations normally found In Rawlnsonde profiles. The smooth- 
ness of the wind profile model is not considered to be a serious deficiency 
because wind biasing Is with respect to the predominant large scale perturba- 
tion in the profile. Thus, even if a single Rawlnsonde profile obtained a 
few hours prior to launch Is used as the basis for wind biasing, the small 
scale perturbations In the profile would be removed before implementation. 
Nevertheless, filtered profiles can still differ greatly from the monthly mean 
profile; therefore, Individual filtered wind profiles that are representative 
of the wind conditions associated with a particular launch would be the most 
desirable basis for wind biasing of launch vehicles. The monthly mean wind 
profile is almost never representative of launch conditions. 

The development of a pre- launch wind monitoring scheme to provide data 
for wind biasing will require knowledge of the change of smoothed wind pro- 
files with -espect to time. This report describes wind bias change with 
respect to time that has been calculated from the VAFB (1965-74) and KSC (1956- 
70) twice daily Rawinsonde series. Each profile In the series was filtered 
before calculation of wind change statistics. The filtering process removed 
the small scale perturbations. Hind change at KSC and VAFB for unfiltered 
profiles has been described In previous reports [1,2]. The methodology used In 
this study is basically the same; wind bias change for time intervals from 0 to 
72 hours at altitudes from 5 to 22 km Is calculated for selected months. Wind 
bias change is presented in terms of statistical summaries of wind bias com- 
ponent change and the modulus of vector wind bias change; the parameters of 
theoretical probability distribution functions representing the wind bias 
change variables are also presented. The validity of the theoretical dis- 
tributions is established by comparing them with the observed distributions. 
These distribution functions can- be utilized to obtain statistical predictions 
of wind change with respect to time. 

This report consists of a Drlef statement of technical background 
(Section II), an analysis of wind change statistics calculated from filtered 
data (Section III) and conclusions (Section IV); the calculated statistics of 
wind bias change with respect to time for selected months 1 km altitude 
Increments from 5 to 22 km to KSC and VAFB are listed in the appendix. 


1 



- . ifcfcMSisa a Ja&>i mzJ &Mi 2J2d jyfii^, ■* frLlyijl ££l 2*jl 'f4 V rtTr^ • I '»)'*} • * 


II. TECHNICAL BACKGROUND 


A. INTRODUCTION 

The large sample of wind profiles obtained at VAFB Is suitable for 
calculation of an equally large sample of wind bias change data. In order to 
readily abstract Information on wind bias change from these data. It Is 
necessary to perform a second series of calculations which provide statis- 
tical summaries of wind bias change. The choice of statistical parameters 
for description of wind bias change Is based in part on the need to specify 
the parameters of theoretical distributions of wind change. These theoretical 
distribution functions are described in detail by Smith (3). The basic 
distribution of the four variables consisting of the zonal and meridional 
components of the wind bias vector at an Initial time and after an elapsed 
time, At, Is quadravarlate normal. The conditional distribution of the wind 
bias components at a specified future time, given the wind bias components 
at an Initial time, Is bivariate normal. The modulus of the wind bias change 
vector is Rayleigh and the distribution of either the zonal or meridional 
wind bias component change Is univariate normal. A significant portion of 
the analytical discussion In Section III of this report Is the presentation 
of observed distributions of wind bias change and comparison with the 
theoretical distributions of wind bias change variables. Succeeding para- 
graphs of this section are concerned with a description of the wind bias 
profile data, the filtering of the data and the definition of statistical 
parameters of wind bias change used in the various theoretical distribution 
functions. 

B. DATA 

The basic winds aloft data are recorded In terms of wind direction, 0 
and magnitude, W. The wind vector Is expressed In the standard meteorological 
coordinate system in which the direction from which the wind Is blowing Is 
measured In degrees clockwise from true north. The zonal component, u, of 
wind vector Is positive for a west (west to east) wind (0-270°) and negative 
for an east (east to west) wind (0-90°) ; the meridional component, v, is 
positive for a south (south to north) wind (0-180° ) and negative for a north 
(north to south) wind (0*0°); u and v are obtained from 0 and W according to: 
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Wind profiles suitable for wind biasing of launch vehicles (defined 
here as wind bias profiles) are calculated by application of an 11-point 
symmetrical Martin-Graham digital low-pass filter to Rawlnsonde profiles. 

The filter removes the small scale perturbations in the wind profile without 
the addition of phase shift to the data. The filter gain and weighting 
functions are listed in Table 1. The effect of the filter on a particular 
profile is illustrated in Figure 1. Application of the filter to a wind pro 
file originally containing data at 1 km intervals from 0 to 27 km produces 
a somewhat abbreviated filtered profile extending from 5 to 22 km. The 
mathematical background and computer code for calculation of the filter gain 
and weighting functions are described by Demandel and Krivo (4). 

The typically large deviation of individual filtered and unfiltered 
Rawlnsonde profiles from the artificial profile composed of the monthly mean 
at each altitude is also illustrated in Figure 1. 

0. DEFINITIONS 


The subscript 0 is used to denote the initial value of a variable, 
and the subscript 1 denotes the variable after an elapsed time. At. Thus: 

Au » Uj - u Q (4) 

Av * Uj - v Q (5) 


Where Au and Av are the components of the wind bias change for a speci- 
fied At. The modulus, R, of the wind bias change with respect to time is given 
by: 

R * J (Au)^ + (Av)^ (6) 
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Table 1. Filter Weights and Effective Response Function of an 
11-Point Martin Graham Low Pass Filter with a Nominal 
Cutoff Frequency of .04 km* 1 and a Termination Frequency 
of .20 km*? 


Filter Weights, h^ 


Response Function, G(f) 


W) 


ML 


L 

.22658165 

.01 

.996 

0 


.02 

.985 

'♦1,-1 

.20033538 

.03 

.966 


.04 

.940 

‘♦2,-2 

. 13609867 

.05 

.905 


.06 

.863 

‘♦3,-3 

.06181594 

.07 

.813 


.08 

.756 

A 

.00628003 

.09 

.693 

tH f -H 


.10 

.624 

‘>5,-5 

.01832084 

.15 

.251 


.20 

-.002 



.30 

.023 



.40 

-.023 



.50 

.023 



.60 

-.023 



.70 

.023 



.80 

-.023 
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Figure 1. Wind Bias Profile Calculated from a Rawinsonde 
Profile and an Artificial Profile Composed of 
Monthly M^ns for the Period 1956-67 at KSC 
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The statistical means are denoted by an overbar, the standard 
deviations and the correlation coefficients are denoted by o x and R(X, Y), 
respectively, with X and Y replaced with the notation appropriate to the 
variable of interest. 

E. STATISTICS 

The wind vector measurements at an Initial time and after an elapsed 
time are treated In this Investigation as a sample from a quadravarlate 
normal distribution defined by the fourteen statistics listed below: 

MEANS 

V V Y 7 1 

STANDARD DEVIATIONS 


CORRELATION COEFFICIENTS 
R(u q , v q ), R(u q , u^) 

R(v 0 . Vj). R(uj, Vj) 

R(Uj, Vq), R(Vj, Uq) 


The fourteen statistics of the quadravarlate normal distribution of 
vector wind difference with respect to time consist of the five bivariate 
normal statistics of vector wind at an Initial time (u n , v ft , o„ o u and 

V U Ua» Vq 

R(uq, Vq)) and the nine statistics Involving component differences which 
can be calculated from the quadravarlate statistics listed above according 
to the following equations: 


MEANS 

Su • ~ u 0 " “l " “o 

_ _ __ ( 8 ) 

AV ■ Vj - Vq ■ Vj - Vq 
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Where R(x,y) Is the correlation coefficient of variables x and y. 
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III. ANALYSIS 


A. INTRODUCTION 

The analysis of wind bias profile change with respect to time follows 
the approach taken In previous studys of vector wind change at KSC and VAFB 
[1,2]. The vectors under consideration have been modified by the filtering 
process described In the previous section. Since the component of wind 
change associated with small scale perturbations In the profile has been 
removed, the calculated wind change Is expected to be smaller for wind bias 
profiles. The objective of this analysis Is the establishment of a theoreti- 
cal basis for estimation of wind bias change. This Is accomplished by 
comparison of theoretical probability distributions which contain wind bias 
change sample statistics as parameters (from the appendix of this report), 
to observed probability distributions of wind bias change. Wind bias change 
with respect to time Is analyzed herein In terms of wind component change, 
unconditional and conditional joint distributions of wind component change, 
and the modulus of vector wind change. 

B. WINO BIAS COMPONENT CHANGE WITH RESPECT TO TIME 

The theoretical probability distribution of wind component change 
with respect to time Is univariate normal with zero mean and standard 
deviation given by Equations 9 and 10; the assumption of zero means of 
component differences is varifled by the sample statistics given In the 
appendix. The theoretical normal distribution of component differences 
can be derived by using either the standard deviations of component differ- 
ences given In the appendix or an estimate which can be obtained from the 
standard deviation of the components If It Is assumed that: 
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Equations 9 and 10 reduce to 


% °U V 1 - «<v V 


CT fc a 
Av 


>1? a v V 1 “ R (V V 


The wind component autocorrelation functions, R(Uj, Uq) and R(v^, Vq) can 
be represented by a negative exponential function of time increment. 


t ; i • e . , 


R(uj. Uq) * EXP ( -bx ) 

R(vj, v 0^ = EXP ^" CT ^ 


where b and c are computed according to 


Zr.ln R^Uj.Uq) 


5 T i ,nR i< v r v o> 


Examples of the decay of the autocorrelation function at 6, 12 and 
18 km during Apr 1 ’! at Cape Kenneviy and Janc'.ry at VAFB are illustrated in 
Figures 2 and 3, respectively; the lines in the figure represent the decay 
rate predicted by Equations 18 and 19. 

Substitution of Equations 18 and 19 into 16 and 17, respectively, 

yields a simple expression for a^ u and o^ y in terms of o u and o y , 
respectively. 


a Au -VT 1 - (-1*7 


° 4V -VT c v V 1 - EXP (-«) 
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Figure 3. Zonal and Meridional Wind Bias Component 

Autocorrelation During January at 6, 12 and 
18 km at Vandenberg AFB (1965-74) 
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Equations 22 and 23 Indicate that and e am asymptotic to 
VF o u and vr o y for large values of t. Thamfom, astlmatas of tha 
extreme value of o^ u and o^ v am obtalnad by setting t aqual to - In 
aquations 22 and 23. 


Tha calculated valuas of b and c at altitudes from 5 to 22 km ovav K3C 
In April and VAFB during January listed In Tabla 2 am also plotted In Figures 
4 and 5. At both locations at altltudas from 5 to 22 km tha decay rate of the 
meridional wind bias component autocorrelation Is larger than the decay rate 
for tr.e tonal component. The variation of the decay rate for the meridional 
component as a function of altitude differs at the two locations. Maximum 
decay for the meridional component occurs at 12 - 13 km over VAftf during Jan- 
uary; In contrast, at kSC during April, the maximum occurs at the extremes 
of altitude range (near 5 and 22 km) and Is a minimum at la km. The decay 
of tonal component autocorrelation decreases steadily with altitude at VAFB 
during January. At KSC during April the decay also decreases with altitude 
but at a very rapid rate and within a restricted altitude range (5 - 15 kml; 
above 15 km the decay Increases. 


The calculated and observed values of o^(r) and o^ y (t) at 6, 12, and 
IS km during April at KSC and January at VAFB are listed In Tables 3 and 4. 

The estimated extreme values of o Uj and o^ y , 

' u VT o y . respective- 
ly), are listed at the bottom of each column of calculated values. The com- 
parisons In Tables 3 and 4 Indicate that o and o^ y can be accurately esti- 
mated by application of Equations 22 and 23, respectively. General application 
of this estimation technique at other locations would require a aw re adequate 
knowledge of the form of the autocorrelation function than Is presently avail- 
able. 


The thvOretlcal distribution of wind bias component differences has 
boon derived from sample estimates of and o^ y and .Yu and ,\v tgiven In 
the appendix 1 for the intervals of 12, 24, 3b and 48 hours at 12 km during 
Ap"i 1 at kSG and January at VAFB; the theoretical normal distributions are 
plotted as straiolit lines In Figures b through 9; the plotted symbols repre- 
sent the observed distributions of Au and ,\v . It is Indicated that the 
ohserv \1 ..istrlbutlon of bias exponent changes is accurately represented 
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Z (km) 

VAFB (JANUARY) 
10 2 b(hr -t ) uAdii -1 ) 

KSC (APRIL) 

10 2 b(hr _1 ) 10 2 c(hr _1 ) 

S 

1.79 

2.72 

1.17 

2.96 

6 

1.80 

2.83 

1.08 

2.64 

7 

1.80 

2.99 

1.04 

2.43 

8 

1.78 

3.16 

1.02 

2.32 

9 

1.72 

3.33 

1.02 

2.25 

10 

1.63 

3.49 

0.99 

2.21 

ii 

1.52 

364 

0.95 

2.15 

12 

1.41 

3.74 

a 89 

2.06 

13 

1 32 

3.77 

0.85 

1.97 

14 

1.25 

3.87 

0.81 

1.89 

IS 

1.21 

3.42 

0.80 

1.86 

16 

1.18 

3.04 

080 

1.89 

17 

1.09 

2.67 

082 

1.99 

18 

0.98 

2.11 

0.84 

2.14 

19 

081 

1.75 

088 

2.38 

20 

070 

1.50 

0.90 

2.79 

21 

065 

1.37 

0.98 

3.71 

22 

083 

1.34 

1.10 

3.34 
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Constants b and c of Equations IS and 19 at Altitudes from 
5 to 2J Km During January at VAFB (1965-’4) and During April 
at Cape Kennedy ( 1956-70) 











Figure 4. Constants b 



t 1 r 

2 3 

10*6 (hr -1 ) 
10* • (hr" 1 ) 


and c of Equations 18 and 19 for Cape 
ng April (1956-70) 
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ALTITUDE 

(KM) 

t 

(HOURS) 



V 

CALCULATED 

OBSERVED 

CALCULATED 

OBSERVED 

18 

12 

3.11 

117 

120 

160 


24 

4.29 

4.13 

126 

189 


38 

113 

104 

4.92 

4.78 


48 

5.78 

5.70 

138 

133 


60 

6.32 

124 

171 

172 


72 

6.78 

159 

195 

199 


oo 

10.04 

— 

— - 

— 

12 

12 

7.00 

123 

7.89 

7.28 


24 

aes 

128 

10.63 

1135 


38 

11.62 

11.50 

1120 

1134 


48 

12.97 

1194 

1136 

1144 


80 

14.16 

14.06 

14.20 

14.23 


72 

15.14 

14.92 

14.)*2 

14.66 


oo 


— 

— — 

— 

8 

12 

6.12 

4.74 

111 

4.76 


24 

7.01 

7.04 

172 

169 


38 

133 

163 

7.69 

7.83 


48 

934 

138 

132 

182 


60 

iai4 

1103 

8.76 

178 


72 

1179 

1154 

106 

182 


oo 

14.68 

— 

— 
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Table 3. 


Calculated [Eqs. 22 and 23] and Observed a . and <ja V 
Bias Profiles During April at Cape Kennedy <lt 6, 12, 


from Wind 
and 18 km 
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ALT1TUDC 

(KM) 

™ "™ " 

r 

(HOURS) 


*\\» 

CALCULATED 

OBSERVED 

CALCULATED 

OBSERV 

It 

12 

3.20 

2,97 

3 56 

2 86 



A 40 

4 08 

4 75 

4 42 


B Kl B 

624 

5,03 

5,49 

5 44 


B 

■-in m 

6 75 

aoi 

6 0 7 


60 

B 

648 

836 

6 43 


72 


7.03 

6 65 

6 65 


A’ 

a 72 

-- 

7 52 


12 

12 

a ie 

8,12 

12 23 

9 9 1 


24 

1243 

12.16 

15,65 

1MU 


36 

1464 

1466 

17 49 

17 56 


46 
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Figure H. Theoretical (Straight Lines) and Observed (Plotted Points) Cumulative Probability 
Distribution of Zonal Wind Bias Component Change, .'.j, with Respect to Time 
Increment, During January at 12 Ion Over Vandenberg AFB (1965-74) 
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by the theoretical normal distribution for a large range of probabilities; 
the deviation of the observed distribution from the theoretical distribution 
at the extreme probabilities Is attributed to the small sample of observations 
and errors In the Rawlnsonde data. 

C. JOINT DISTRIBUTION OF WIND BIAS COMPONENT CHANGES WITH RESPECT TO TIME 

The joint distribution of zcnal and meridional wind bias component 
change with respect to time (Au and Av) can be approximated by a bivariate 
normal distribution. A useful property of such a distribution Is that an 
ellipse can be calculated which contains the end points of a specified per- 
cent of vectors having components Au and Av. A detailed description of the 
derivation of probability ellipses and plotting methodology is given by Smith 
(3). The five parameters of the bivariate normal distribution of Au and Av, 
calculated for each monthly reference period at KSC and VAFB at 1 km altitude 
intervals from 5 to 22 km are listed In the appendix. 

The degree of approximation of the bivariate normal distribution to 
the observed distribution can be evaluated by comparison of the observed 
percentage of vectors which are contained within the ellipse to that pre- 
dicted by the ellipse at a specified probability level. For example, for 
a sample of 1,000 vectors, 950 of the vectors should terminate within the 
95 percent (theoretical P ■ .95) ellipse calculated from the bivariate 
statistics of the 1,0C0 vectors; however, a plot of the 1,000 vectors 
could Indicate that only 945 vectors (observed p045) terminate within the 
95 percent ellipse. For illustration on a linear graph comparison of the 
theoretical to the observed P is given in terms of the parameter X given 
by 

x e - Vf V- In (1-P) ( 24 ) 

A comparison of theoretical and observed values of X„ at 12 km during 
April at KSC and January at VAFB for time Intervals of 12, 24, 36 and 48 hours 
is Illustrated In Figures 10 and 11. Perfect agreement between theoretical 
and observed X g Is represented by a line drawn from the origin with a slope, 

B, equal to 1. The calculated least squares slopes are given in the figure 
legend. The plots indicate a tendency for the theoretical X g tc exceed the 
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Figure 10. Observed as a Function of Theoretical for a HI variate 
Normal Distribution of Wind Bias Component Changes (An, AvI 
with Respect to Time at It km During April at 

Cape Kennedy 
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observed for large values of X^. The Interpretation of this result Is 
that for extreme probabilities the theoretical distributions predict fewer 
wind change vectors terminating outside the ellipse than Is observed. These 
results may have to be taken Into consideration if engineering application of 
theoretical wind change statistics beyond 95 percent level Is required. 

The 95 percent probability ellipses for the joint distribution of 
wind bias component changes with respect to time at 6, 12, and IS km during 
April at KSC and January at VAFB are Illustrated in Figure 12; the relatively 
large changes with respect to time at 12 km Is clearly Illustrated at both 
locations. 

0. MODULUS OF VECTOR WIND BIAS CHANGE WITH RESPECT TO TIME 

If wind bias change with respect to time has a distribution which Is 
bivariate normal, the modulus R, of the wind bias change vector (defined by 
Equation 5) has a Rayleigh distribution. Since the Rayleigh distribution 
cannot be Integrated In closed form, numerical integration 1 c required to 
obtain the cumulative probability distribution. Derivation of the Rayleigh 
distribution, given the five bivariate normal distribution statistics, re- 
quires summation involving products of the modified Bessel function of the 
first kind. Smith i 31 summarizes the basic equations for the Rayleigh 
distribution derived by Wler (51 and extended by Yadavalll (61 to include 
the condition for correlated variables. The Rayleigh distribution reduces 
to the Integrable classical form If it is assumed that the components of 
the vector wind change are Independent and that they have zero means and 
equal standard deviations; the classical Rayleigh probability density 
function is 


f(R) - A EXP (-RV2 a 2 ) R • 0 (25) 

o* - 

Integration of Equation 25 from zero to a specified value of R yields 
the cumulative probability that R v R* where, 

Pr (R < R*) » 1 -EXP (-RV2o 2 ) R • 0 (26) 

where o ■ o. ■ o. 

Au Av 
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Figure 12. Joint Distribution of 95 Percent Wind Bias 

Component Changes with Respect to Time at 6, 12 
and 18 km During April at Cape Kennedy (1956-70) 
and January at Vandenberg (1965-74) 
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Since the standard deviation of the component difference can be 
expressed as a function of the standard deviation of the components 
(Equations 22 and 23) it follows that 


Pr [R < R*J - 1 - EXP 



R- 

(1-EXP (-kr)i 


( 27 ) 


where o k and k correspond to either o u and b or o y and C given in Equations 
22 and 23. 


An expression for R given a particular probability, Pr [R ^ R*i , is 
obtained by solution of Equation 27 to obtain 

R * ^ Vk Vi - EXP (-kt) CCS) 

where Is derived from Equation 24 denoting Pr (R < R*] by P 

The choice of o k » o y and k « c (from Equation 23) at 12 km during 
April at (CSC and January at VAFB yields the most accurate approximation of 
the cumulative Rayleigh distribution obtained by numerical integration of 
Equation 28 in Reference 1. Comparisons of the 99, 95, and 50 percentile 
modulus of the wind change vector with respect to time based on the Rayleigh 
(Equation 28, Reference 1) and the classical Rayleigh (Equation 27) are 
Illustrated in Figures 13 and 14; the rather good agreement between the 
distributions for time intervals from 12 to 72 hours is attributable to the 
accuracy of the simplifying assumptions described above. There is a slight 
tendency, especially for time intervals 06 hours for the classical Rayleigh 
to be larger than the Rayleigh. 

The remaining question is: How well do these theoretical distributions 
compare with observed observations? Comparisons of observed and theoretical 
values of R for time intervals of 12, 24 36 and 48 hours at 12 km during 
November, December and January at VAFB are given in Tables 5 and 6; column 
II of the tables contains R calculated according to the classical Rayleign 
distribution with o equal to the monthly value of o at 12 km and k 
equat to the decay constant, c, in the monthly exponential least squares 
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fit to the v component autocorrelation function (Equation 23); Column I was 
obtained by numerical Integration of the Rayleigh distribution. It Is 
Indicated that the observed cumulative distribution agrees fairly well with 
the theoretical distribution for probabilities less than .95 to .975. For 
large probabilities, there is a consistent tendency for the theoretical dis- 
tribution to underestimate the observed distribution. This tendency is 
attributable to the small sample of data avilable at the extreme probabil- 
ities and errors in the Rawinsonde data. 

E. CONDITIONAL VECTOR WIND BIAS ELLIPSES 

Prior knowledge that environmental constraints necessary to assure 
the success of a space vehicle launch will be satisfied implies that there 
is a capability for prediction of environmental parameters; the prediction 
can be based on knowledge of conditions prior to launch. With regard to 
winds aloft, prior conditions are typically based on Rawinsonde or Jimsphere 
wind profiles. A typical question that could be posed before launch Is: 
Given a measurement of the wind bias vector 12 prior to launch at 12 km, 
will the wind bias vector at launch time be within 95 percent reference 
month wind ellipse? A question of this type can be answered If the distri- 
bution of vector wind bias components at an Initial time, T Q , and at a 
future time, Tj, can be approximated by a quadravariate normal distribution. 
Given the components of the bias vector at Tq, the conditional distribution 
of the bias vector at T^ is bivariate normal. Smith tl] describes the deri- 
vation of the conditional bivariate normal distribution and documents the 
computer program used in this investigation for calculation of these dis- 
tributions. Figures 15 and 16 illustrate the 95 percent conditional bi- 
variate normal distributions at 12 km that have been calculated for time 
Increments of 12, 24, 36, 48, 60 and 72 hours for the month of April at KSC 
and January at VAFB . Five vectors were selected as given Initial conditions 
for calculations of the conditional ellipses. The components of the vectors 
are defined below: 

1. Monthly bias component means 

2. Maximum zonal wind bias and the corresponding meridional wind 
bias from the monthly 95 percent vector wind bias ellipse. 

3. Minimum zonal wind bias and the corresponding meridional wind 
bias from the monthly 95 percent vector wind bias ellipse. 
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Figure 15. April Conditional 95 Percent Wind Bias Ellipses at 12 km 
for Time Increments of 12, 24, 36, 48, 60 and 72 Hours at 
Cape Kennedy (1956-70) 










4 . 


Maximum meridional wind bias and the corresponding zonal 
wind bias from the monthly 95 percent vector wind bias 
ellipse. 

5. Minimum meridional wind bias and the corresponding zonal wind 
bias from the monthly 95 percent vector wind bias ellipse. 

The given vectors are specified In the inset of Figures 15 and 16 (polar 
form, at 12 km) and In Tables 7 and 8 (component form, at 6, 12, 18 km). 

The conditional ellipses Illustrated at the center of Figures 15 
and 16 show that If the observed wind vector has components equivalent to 
the monthly mean bias components (condition 1) then 95 percent of the wind 
vectors after elapsed times as large as 72 hours will fall within the monthly 
95 percent ellipse. Therefore, satisfaction of a launch constraint which 
states that the wind bias vector must be included within *he 95 percent 
monthly ellipse would be assured for periods as long as 72 hours following 
an observation of a wind vector having components which correspond to the 
monthly means. The conditional ellipses based on selection of given wind 
bias vectors that terminate on the monthly 95 percent ellipse (conditions 
2 through 5) have a significant proportion of their area lying outside the 
monthly 95 percent ellipse; as the time increment increases this proportion 
decreases but remains significant for a time increment as large as 72 hours. 
This implies that a significant proportion of wind bias vectors will not 
satisfy a launch constraint oased on the 95 percent wind bias ellipse for 
periods as long as 72 hours (or longer if these calculations are extended) 
following an observation of a wind bias vector which terminates on the 95 
percent ellipse. 

The wind direction characteristics of a wind bias ellipse can be 
described in terms of the angles associated with wind bias vectors constructed 
between the origin and the center ol ‘.he ellipse (at the component means) and 
between the origin and the two tangent positions to the ellipse. The three 
vectors constructed in this manner and the angles t? A , c?g, and At? are 
illustrated in Figure 17, the range of wind angles, e R , is t? A to t?g. The 
angles 6 R ,t? E , At? calculated from five 95 percent conditional ellipses for 
April at KSC and January at VAFB at 6, 12 and 18 km are listed in Tables 7 
and 8. 
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CONDITION (*) 

0. 

°A 

*B 

M 


1 

278 

248 

313 

65 


2 

2b6 

262 

281 

39 


3 

20 

• 

• 
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4 

244 

225 

266 

41 


5 

323 

293 

349 

56 



Cl THE FIVE CONDITIONAL DISTRIBUTIONS ARE BASED ON THE FIVE GIVEN WINO VECTORS 
LISTED BELOW. CONDITION 1 IS BASED ON MONTHLY MEAN WINO BIAS COMPONENTS FOR 
THE PERIOD 1956 -70; CONDITIONS 2 THRU 5 ARE FROM THE 95 PERCENT VECTOR WIND 
BIAS ELLIPSE AT 12 KM THAT WAS CALCULATED FROM TWICE DAILY FILTERED RAWINSONDE 
OATA DURING THE PERIOD 1966-70. 

* 95 PERCENT CONDITIONAL ELLIPSE COVERS ALL QUADRANTS. 


CONDITIONS 



M/SEC 

1 

U, V 

30.44 

- 4.35 

2 

u mw * 

68.53 

4.46 

3 

u mln* ¥ 

- 7.66 

- 13.18 

4 

“• *ro«x 

41.94 

24.83 

5 

u > 

18.94 

- 33.53 
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Table 7. Wind Direction (Decrees) Characteristics of 95 Percent 

Conditional Vector Wind Bias Ellipses at i2 km Over Cape 
Kennedy During April for an Elapsed Time, x, of 12 Hours 
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THE FIVE CONDITIONAL DISTRIBUTIONS ARE BASED ON THE FIVE GIVEN WIND VECTORS 
LISTED BELOW. CONDITION 1 IS BASED ON MONTHLY MEAN WIND BIAS COMPONENTS FOR 
THE PERIOD 1966-74; CONDITIONS 2 THRU 5 ARE FROM THE 96 PERCENT VECTOR WIND 
BIAS ELLIPSES AT 12 KM THAT WERE CALCULATED FROM TWICE DAILY FILTERED 
RAWINSONDE DATA OURING THE PERIOD 1966 -74 
98 PERCENT CONDITIONAL ELLIPSE COVERS ALL OIJADRANTS. 


CONDITION 


M/SEC 

1 

u, V 

2138 

- 495 

2 

U V 

max' w 

63.60 

4.47 

3 

u V 

min* 

- 18.78 

- 14.37 

4 

u - v m« 

34.11 

30.26 

5 

u - *min 

12.61 

- 4015 
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Table 8. Wind Direction (Degrees! Characteristics of os Porcent 

Conditional Vector Wind Bias ellipses at 1C km Over Vandonberg 
AFB During January for an t lapsed Time, i, of 1C Hours 
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IV. CONCLUSIONS 


The analysis presented In the preceding section for selected months 
and altitudes illustrates how various theoretical distribution functions 
can be used for calculation of wind bias change with respect to time at 
Cape Kennedy, Florida and Vandenberg AFB, California. The calculations can 
be made by utilization of the statistics given in the appendix for any refer- 
ence month at 1 km altitude increments from 0 to 27 km. It also has been 
shown that the techniques originally used to describe wind change observed 
in unfiltered Rawlnsonde profiles can also be applied with equivalent accuracy 
to describe wind bias change. 

The basic underlying assumption for the calculation of the distribu- 
tions of wind bias change is that the joint distribution of the four variables 
represented by the components of the wind bias vector at any initial time 
and after a SDecified elapsed time is quadravariate normal. If the wind bias 
vector is specified at an initial time, then the conditional joint distri- 
bution of the wind bias components at a future time is bivariate normal. 

Since each of the variables of the quadravariate normal distribution is 
normal and the difference of two normal distributions is normal, it follows 
that wind bias component change is also normal and the joint distribution 
of zonal and meridional wind bias change is bivariate normal. The modulus 
of bivariate normally distributed variables has a Rayleigh distribution. 
Therefore, the modulus of vector wind bias change with respect to time is 
Rayleigh. 

Sample distributions based on reference month Rawinsonde data obtained 
during January 1965-74 at Vandenberg AFB and April 1956-70 at Cape Kennedy 
agree reasonably wel 1 with the aforementioned theoretical distributions. 

The standard deviation of wind bias component change with respect to 
time is the only statistic required for determi nation of the theoretical 
probability distribution (normal with zero mean) of wind bias component 
change. It has been shown than over a large range of altitudes that this 
statistic can be estimated from wind bias component standard deviation and 
the decay constant of the component theoretical autocorrelation function 
(Table 1). The assumption of exponential decay of the autocorrelation 
function is reasonably accurate in most instances to time Increments as large 
as 72 hours at both locations. 


* u— 


v 



T'l I Lt I 1,1:1 


The observed modulus of vector wind bias change with respect to time 
Is systematically larger than the predicted modulus (Section III.C) for 
extreme probabilities. This may be attributable to Inadequacy of the theory 
or Inaccuracies of the data which affect the observed distribution at the 
extreme probabilities. If the theoretical distribution at extreme proba- 
bilities Is to be used In engineering applications, it will be necessary to 
explain these systemmatic differences. 
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APPENDIX 

This appendix contains two sets of reference month quadravarlate and 
conditional bivariate normal statistics of variables X, Y, XP and YP, at 
1 km Intervals from 5 to 22 km. The statistics were calculated from serially 
complete twice dally wind bias profiles calculated from Rawlnsonde profiles 
obtained during the period 1965-74 at VAFB and 1956-70 at KSC. The notation 
for the variable given In Section II of this report differs from the notation 
established for the computer output given herein; the notations are compared 
In Table A-l. 


TABLE A-l. NOTATION OF VARIABLES 


Computation Set 

A 

B 

Variable 

Text 

(Sect. II) 

Computer 

Output 

Text 

(Sect. II) 

Computer 

Output 

X 

u 0 

u(at T) 

u 0 

u(at T) 

Y 

v 0 

v(at T) 

v 0 

v(at T) 

XP 

U 1 

u(at T+DT) 

U 1 • u 0 

» Au 

u(at T+DT) 
-u(at T) 

YP 

V 1 

v(at T+DT) 

V 1 ‘ v 0 
* Av 

v(at T+DT) 
-v(at T) 


Table A-l shows that the quadravarlate statistics of computation set "A" 
are for wind bias components at an Initial time and after a specified time 
increment; the statistics for set "B" are for wind bias components at an 
initial time and wind bias component change after a specified time increment. 
The reference month quadravarlate normal statistics at a particular altitude 
for six time increments (12, 24, 36, 48, 60 and 72 hours) are listed in the 
lower left of each page of computer listing; the six sets of conditonal 
bivariate normal statistics corresponding to the six time Increments are 
listed in the lower right. The data were conditioned on monthly means 
for the entire data sample. The derivation of the conditional bivariate 
normal statistics for any other given vector involves recalculation of the 
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conditional means according to equations A-l and A-2; the standard de- 
viations and correlation coefficients do not have to be recalculated 
because they are independent of the given wind vector. 


x c jxp*«x + 


[ (R(x,xp) - R(x,yp) R(xp ,yp) ) (xp* - xp) (o x /o xp ) 
♦ (R(x,yp) - R(x,xp) R(xp,yp) ) (yp* - yp) {c^/c )l 

— — 2 

1 - [R(xp,yp)j 


(A-l) 


[(R(y.xp) - R(y,yp) R(xp.yp)) (xp* - xp) (■? /o ) 

_ _ + (R(y.yp) - R(y.xp) R(xp.yp)) (yp* - yp' (o /o )] 

ylyp*-y * 1 — ^ — 

c 1 - [Rixp.ypir 


( A-2 ) 


where, x c and y c are the mean components of the conditional distribution, 

xp* and yp* are the components of the given vector and 

o ,o ,o and o are equivalent to S.D.x, S.D.y, S.D.xp and 
x y xp yp 

S.D.yp, respectively given in the computer listings. 
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QUADRA VAR IMF NORMAL STATISTICS OF *,V,*P,VT> • CONDITIONAL BIVARIATE NORMAL STATISTICS 
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